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(54) FLOW RATE CALCULATING DEVICE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a flow rate calculating 
device capable of precisely calculating the flow rate of fluid 
passing through a ventilation path. 

SOLUTION: This flow rate calculating device calculates the air 
flow rate passing through a throttle valve 23 provided on the 
way to an intake pipe 20. On the assumption that the flow rate 
of fluid passing the throttle valve 23 is set to mt, the opening 
area of the throttle valve 23 is set to Ad, an upstream pressure 
in the upstream side of the throttle valve 23 is set to Pu, an 
upstream density in the upstream side of the throttle valve 23 
is set to pu, a downstream pressure in the downstream side of 
the throttle valve is set to Pd, and a specific heat ratio of the 
passing air is set to k, the air flow rate mt is calculated from 
the following expression. mt= Ad.(Pu.pu)1/2.((((k-1)/(2.k)).(1- 
Pd/Pu)+Pd/Pu).(1 -Pd/Pu))1 /2. 
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* NOTICES * 

iTPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 

[Claim 1] In the flow rate calculation equipment which computes the fluid flow which passes the adjustable 
converging section formed in the middle of an aeration path The opening area of mt and said adjustable 
converging section for said fluid flow which passes said adjustable converging section Ad, The upper 
consistency of the upstream of said adjustable converging section [ in / for the upper pressure of the 
upstream of said adjustable converging section in said aeration path / Pu and said aeration path ] rhou, Flow 
rate calculation equipment characterized by having an operation means to compute said fluid flow mt by the 
following formulas when the ratio of specific beat of Pd and said fluid is set to k for the down-stream 
pressure of the downstream of said adjustable converging section in said aeration path. 
mt=Ad (Pu-rhou), 1/2 (k-1) (/(2andk)) (-(1-Pd/Pu)+Pd/Pu) (- (1-Pd/Pu)), 1/2. 

[Claim 2] Said operation means is flow rate calculation equipment according to claim 1 characterized by 
computing said fluid flow mt using the following formulas. 
mt=Adand(Pu-rhou) 1/2, andphi (Pd/Pu). 

In addition, Function phi (Pd/Pu) is 1/2 at the time (k-1) of 1/2 and (Pd/Pu)> (l/(l+k)) (/(2andk)) (-(1- 
Pd/Pu)+Pd/Pu) (- (1-Pd/Pu)) at the time (k/(2- (k+1))) of (Pd/Pu)<= (l/(l+k)). 

[Claim 3] Said adjustable converging section is flow rate calculation equipment according to claim 1 or 2 
characterized by being the throttle valve prepared in an internal combustion engine's inhalation-of-air path. 
[Claim 4] Said adjustable converging section is flow rate calculation equipment according to claim 1 or 2 
which is a duty control valve and is characterized by setting up said opening area based on duty ratio. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the flow rate calculation equipment which computes the 
fluid flow which passes the adjustable converging section formed in the middle of an aeration path. 
[0002] 

[Description of the Prior Art] It is the inhalation air content calculation approach which computes the air 
content inhaled by an internal combustion engine's cylinder through the inlet pipe equipped with the throttle 
valve, and it considers that a throttle valve is an orifice and what calculates the air content which passes 
through that according to a throttle opening using the formula of hydrodynamics is known so that JP,6- 
74076,A may indicate an aeration path conventionally as the technical technique which computes the flow 
rate of flowing fluid and which carries out flow rate calculation. 

[0003] This flow rate calculation approach presupposes that it computes so that next presumed calculation 
may not be influenced, even if there is that presumed error while raising the presumed precision of a flow 
rate using a real pressure on the occasion of flow rate calculation. 
[0004] 

[Problem(s) to be Solved by the Invention] however, the fluid which uses the relational expression of 
adiabatic change and circulates by this flow rate calculation approach — setting — etc. — from [ computing 
the flow rate on the assumption that entropy conditions are satisfied ] — etc. — there is a trouble that an exact 
flow rate is uncomputable under the situation that entropy conditions are not satisfied. Since ** entropy 
conditions are not necessarily satisfied when computing the air content which flows an internal combustion 
engine's inlet pipe especially, calculation of an exact flow rate is difficult. 

[0005] Then, this invention is made in order to solve such a trouble, and it aims at offering the flow rate 
calculation equipment which can compute the flow rate of flowing fluid for an aeration path to accuracy. 
[0006] 

[Means for Solving the Problem] Namely, the flow rate calculation equipment concerning this invention is 
set to the flow rate calculation equipment which computes the fluid flow which passes the adjustable 
converging section formed in the middle of an aeration path. The opening area of mt and an adjustable 
converging section for the fluid flow which passes an adjustable converging section Ad, When the ratio of 
specific beat of Pd and a fluid is set to k, the down-stream pressure of the downstream of an adjustable 
converging section [ in / for the upper consistency of the upstream of an adjustable converging section / in / 
for the upper pressure of the upstream of the adjustable converging section in an aeration path / Pu and an 
aeration path / rhou and an aeration path ] It is characterized by having an operation means to compute the 
fluid flow mt by the following formulas. mt=Ad (Pu-rhou), 1/2 (k-1) (/(2andk)) (-(1-Pd/Pu)+Pd/Pu) (- (1- 
Pd/Pu)), 1/2. 

[0007] Moreover, the flow rate calculation equipment concerning this invention is characterized by the 
above-mentioned operation means computing the fluid flow mt using the following formulas. mt=Adand 
(Pu-rhou) 1/2, andphi (Pd/Pu). In addition, Function phi (Pd/Pu) is 1/2 at the time (k/(2~ (k+1))) of (Pd/Pu) 
<= (l/(l+k)), and is 1/2 at the time (k-1) of (Pd/Pu)> (l/(l+k)) (/(2andk)) (-(1-Pd/Pu)+Pd/Pu) (- (1-Pd/Pu)). 
[0008] Moreover, the flow rate calculation equipment concerning this invention is characterized by the 
above-mentioned adjustable converging section being the throttle valve prepared in an internal combustion 
engine's inhalation-of-air path. 

[0009] Moreover, the above-mentioned adjustable converging section is a duty control valve, and the flow 
rate calculation equipment concerning this invention is characterized by setting up opening area based on 
duty ratio. 
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[0010] According to these invention, since the fluid flow is computed without being based on the relational 
expression of adiabatic change, calculation of the fluid flow can compute to accuracy under the conditions 
on which entropy conditions are not satisfied. 
[0011] 

[Embodiment of the Invention] Hereafter, the gestalt of operation of this invention is explained to a detail 
with reference to an accompanying drawing. In addition, in explanation of a drawing, the same sign is given 
to the same element, and the overlapping explanation is omitted. 
(The first operation gestalt) 

[0012] The explanatory view of the flow rate calculation equipment applied to this operation gestalt at 
drawing 1 is shown. 

[0013] As shown in this Fig., the flow rate calculation equipment concerning this operation gestalt is applied 
to the calculation of an air content which passes the throttle valve 23 of the inlet pipe 20 of an engine 2. The 
engine 2 is equipped with the adjustable valve timing device 5 in which have an adjustable valve gear, for 
example, the closing motion timing of an inlet valve 3 and an exhaust valve 4 is changed. 
[0014] It connects with ECU6 electrically, and the adjustable valve timing device 5 operates based on the 
control signal outputted from ECU6, and outputs the detecting signal about valve timing to ECU6 through 
the detection sensors 7, such as a cam position sensor. 

[0015] The crank position sensor 12 is formed in the engine 2. The crank position sensor 12 is a sensor 
which detects an engine speed, and it connects with ECU6 and it outputs a pair necropsy appearance signal 
to ECU6. 

[0016] The injector 9 which injects a fuel is formed in the combustion chamber 8 at the engine 2. An 
injector 9 is a fuel-injection means to supply a fuel to a combustion chamber 8, and is installed every 
cylinder 10 with which an engine 2 is equipped. The combustion chamber 8 is formed above the piston 1 1 
arranged in the cylinder 1 0. The inlet valve 3 and the exhaust valve 4 are arranged in the upper part of a 
combustion chamber 8. 

[0017] The inlet pipe 20 which consists of an intake manifold etc. is connected to the upstream of an inlet 
valve 3. The air cleaner 22 is installed in the upstream of an inlet pipe 20. Moreover, the throttle valve 23 is 
formed in the middle of the inlet pipe 20. 

[0018] A throttle valve 23 is an adjustable converging section formed in the middle of an inhalation-of-air 
path, operates based on the control signal of ECU6, changes a throttle opening, and changes the opening 
area of an inhalation-of-air path in connection with it. The throttle opening of a throttle valve 23 is detected 
by the throttle position sensor 24, and is inputted into ECU6. 

[0019] The air flow meter 25 is formed in the down-stream location of an air cleaner 22. An air flow meter 
25 is an inhalation air content detection means to detect an inhalation air content. The detecting signal of an 
air flow meter 25 is inputted into ECU6. 

[0020] ECU6 controls the whole equipment of flow rate calculation equipment, and the computer containing 
CPU, ROM, and RAM is constituted as a subject. The various control routines which contain a flow rate 
calculation routine in ROM are memorized. 

[0021] Next, actuation of the flow rate calculation equipment concerning this operation gestalt is explained. 

[0022] Drawing 2 is a flow chart which shows actuation of flow rate calculation equipment. 

[0023] Step S10 (it is only hereafter indicated as "S10".) of this Fig. Suppose that it is the same about other 

steps. Reading of the throttle opening TA, engine-speed NE, and valve timing VT is performed. 

[0024] Reading of the throttle opening TA is performed based on the output signal of a throttle position 

sensor 24. Reading of an engine speed NE is performed based on the output signal of the crank position 

sensor 12. Reading of valve timing VT is performed based on the output signal of the detection sensor 7. 

[0025] And it shifts to SI 2 and the operation of the air flow rate mt which passes the throttle valve 23 of an 

inlet pipe 20 is performed. The operation of this air flow rate mt is performed based on the following 

formula (1). 

[0026] 

mt=Adand(Pu-rhou) 1/2, andphi (Pd/Pu) .... (1) 

[0027] Here, Function phi (Pd/Pu) is 1/2 at the time (k/(2- (k+1))) of (Pd/Pu)<- (l/(l+k)), and are the time 
(k-1) of (Pd/Pu)> (l/(l+k)) (/(2andk)) (-(1-Pd/Pu)+Pd/Pu) (1-Pd/Pu), and 1/2. 

[0028] Moreover, it is the ratio of specific beat of the air with which the upper pressure of the upstream of 
the throttle valve 23 of an inlet pipe 20 and rhou are inhaled for opening area [ in / in Ad / the location of the 
throttle valve 23 of an inlet pipe 20 ], and Pu, and the down-stream pressure of the downstream of the 
throttle valve 23 of an inlet pipe 20 and k are inhaled for the upper consistency of the upstream of the 
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throttle valve 23 of an inlet pipe 20, and Pd. 

[0029] The operation of the concrete air flow rate mt is performed as follows. 

[0030] Since it is decided that an air flow rate and the pressure-of-induction-pipe force will be meaning 
when an engine 2 is a steady state, if the amount of steady flows at this time is set to mtTA and static 
pressure is set to PdTA, the following formula (2) will be materialized based on an above-mentioned 
formula (1). 
[0031] 

mtTA=Adand(Pu-rhou) 1/2, andphi (PdTA/Pu) .... (2) 

[0032] Here, if it supposes that a formula (1) is materialized also in transients other than a steady state, and 
is dealt with as that for which Ad (Pu-rhou) and 1/2 do not depend on up-and-down differential pressure by 
the same throttle opening in a formula (1) and a formula (2) and Ad (Pu-rhou) and 1/2 are eliminated by the 
formula (1) and the formula (2), the following formula (3) will be obtained. 
[0033] 

Mt=mtTA-phi(Pd/Pu)/phi (PdTA/Pu) .... (3) 

[0034] And since the throttle opening TA, engine-speed NE, and valve timing VT can determine mtTA at 
the time of a steady state, and phi (PdTA/Pu), mtTA/phi (PdTA/Pu) is set up as a table of the throttle 
opening TA, engine-speed NE, and valve timing VT, and phi (Pd/Pu) is set up as a table of Pd. 
[0035] Thereby, according to the throttle opening TA, engine-speed NE, and valve timing VT, the air flow 
rate mt which passes a throttle valve 23 is computable using a formula (3). 

[0036] In addition, in calculation of the above-mentioned air flow rate mt, the air flow rate mt of a transient 
may be computed by setting up mtTA and PdTA at the time of a steady state as a table of the throttle 
opening TA, engine-speed NE, and valve timing VT, respectively, setting up phi (Pd/Pu) as a table of Pd, 
and referring to by Pd and PdTA. 

[0037] Moreover, since the inhalation air content to the throttle opening TA and an engine 2 is equal in the 
case of a steady state, the air flow rate mt of a transient may be computed by setting up PdTA as a table of 
the throttle opening TA, engine-speed NE, and valve timing VT, computing mtTA, setting up phi (Pd/Pu) as 
a table of Pd, and referring to by Pd and PdTA with a charging efficiency. 

[0038] Moreover, the upper pressure Pd may form a pressure sensor in an inlet pipe 20, may use an actual 

measurement for it, and may determine it by presumption from detection values, such as an air flow meter 

25. Furthermore, the table of phi (Pd/Pu) may set an argument to Pd, or may make it Pd/Pu. 

[0039] Next, the calculation result in the flow rate calculation equipment concerning this operation gestalt is 

explained. 

[0040] The calculation result of the throttle passage air content in the flow rate calculation equipment 
applied to this operation gestalt at drawing 3 , the calculation result of the throttle passage air content in the 
flow rate calculation equipment used as the example of a comparison, and the location survey result of a 
high response flowmeter are shown. Time amount change of a throttle opening [ in / in drawing 4 / 
calculation of drawing 3 ] and drawing 5 express time amount change which is the pressure-of-induction- 
pipe force (intake pressure). 

[0041] As shown in drawing 3 , it turns out that the calculation result in the flow rate calculation equipment 
concerning this operation gestalt is well followed in the actual measurement of a high response flowmeter, 
and the flow rate is computed by accuracy. On the other hand, it turns out that the calculation result in the 
flow rate calculation equipment used as the example of a comparison is not followed in the part into which 
the passage air content started at the actual measurement of a high response flowmeter, and the flow rate is 
not computed by accuracy. 

[0042] Here, what computes the throttle passage air content mt based on the following formula (4) was used 

for the flow rate calculation equipment used as the example of a comparison. 

[0043] 

mt=Ad (Pu-rhou), 1/2, phi 1 (Pd/Pu) .... (4) 

[0044] Functions phi 1 (Pd/Pu) are 1 (2/(k+l))/(k-l) (2.k/(k+l)) and 1/2 here at the time of (Pd/Pu)<=(2/ 
(k+1)) k/(k-l). It is 1 (- (2 / k-(Pd/Pu) (k+l)/k) (2.k/(k-l)) (Pd/Pu))/2 at the time of (Pd/Pu)>(2/(k+l)) k/(k- 

[0045] the formula (4) of calculation of this example of a comparison — a law of mass conservation (formula 
(5)) and the law of conservation of energy (formula (6)) ~ reaching — etc. — it is the formula drawn by the 
relational expression (formula (7)) of the adiabatic change on condition of entropy condition formation. 
[0046] 

Au-rho u-vu=Ad-rho d-vd .... (5) 
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[0047] 

vu2/2+(k/(k-l)) -(Pu/rhou) =vd2/2+(k/(k-l)) - (Pd/rhod) .... (6) 
[0048] Pu/rho uk=Pd/rho dk .... (7) 

[0049] In addition, rhod is [ the air rate of flow of the upstream of a throttle valve 23 and vd of the down- 
stream consistency of the downstream of the throttle valve 23 of an inlet pipe 20 and vu ] the air rates of 
flow of the downstream of a throttle valve 23. 

[0050] On the other hand, the flow rate calculation equipment concerning this operation gestalt is computing 
the throttle passage air content mt based on an above-mentioned formula (1). A formula (1) is a formula 
drawn by an above-mentioned law of mass conservation (formula (5)), the law of conservation of energy 
(formula (6)), and the following law of conservation of momentum (formula (8)). 
[0051] 

rhod-vd2, Ad-rhou-vu2, andAu=Pu-Au-Pd-Ad+p- (Ad-Au) .... (8) 

[0052] In addition, p is the whole upstream mean pressure of a throttle valve 23. 

[0053] If the opening area Au of the upstream is considered as infinity when it considers as the opening area 
Au of the upstream of a throttle valve 23, consistency rhou, a pressure Pu, and the rate of flow vu and 
considers as the opening area Ad of a throttle valve 23, consistency rhod, a pressure Pd, and the rate of flow 
vd, as shown in drawing 6 , it will be set to vu=0 from the law of mass conservation of a formula (5). 
[0054] Moreover, if vu=0 is taken into consideration, it will become p=Pu from the law of conservation of 
momentum of a formula (8). Therefore, a law of conservation of momentum serves as rhod-vd2=Pu-Pd. 
[0055] Therefore, a formula (1) and (the case of (Pd/Pu)> (l/(l+k))) are drawn as a throttle passage air 
content mt. Drawing 7 graph-izes phi (Pd/Pu) of the formula (1) in the flow rate calculation equipment 
concerning this operation gestalt (in the case of k= 1 .4). 

[0056] As mentioned above, according to the flow rate calculation equipment concerning this operation 
gestalt, the air flow rate which passes the throttle valve 23 of an inlet pipe 20 is computed using the formula 
(1) drawn by a law of mass conservation (formula (5)), the law of conservation of energy (formula (6)), and 
the law of conservation of momentum (formula (8)), without being premised on ** entropy condition 
formation. For this reason, the throttle- valve 23 passage air flow rate in which ** entropy conditions are not 
necessarily satisfied is computable to accuracy. 
(The second operation gestalt) 

[0057] Next, the flow rate calculation equipment concerning the second operation gestalt is explained. 
[0058] Although the flow rate calculation equipment concerning the first operation gestalt explained the 
case where it applied to the calculation of an air content which passes the throttle valve 23 of the inlet pipe 
20 of an engine 2, the flow rate calculation equipment concerning this operation gestalt is flow rate 
calculation equipment which computes the fluid flow which passes the adjustable converging section 
formed in the middle of an aeration path, and when an adjustable converging section is a duty control valve, 
it is applied. 

[0059] For example, control valves, such as an EGR valve and a purge control valve, repeat ON and OFF, 
and the flow rate calculation equipment concerning this operation gestalt is applied to what computes the 
flow rate which circulates using the formula (1) mentioned above when duty control was carried out and the 
fluid flow was adjusted. 

[0060] In this case, when the control frequency of duty control is smaller than the count frequency of a flow 
rate, an opening area actual as an opening area Ad of a formula (1) is substituted, and calculation of a flow 
rate is performed. 

[0061] On the other hand, as shown in drawing 8 , when the control frequency of duty control is larger than 
the count frequency of a flow rate, the opening area Ad of a formula (1) is set up based on the duty ratio of 
duty control. That is, the average opening area in computation time spacing deltat is used as an opening area 
Ad, and calculation of a flow rate is performed. 

[0062] Even if it is in the flow rate calculation equipment concerning such this operation gestalt, calculation 
of a flow rate is possible like the flow rate calculation equipment concerning the first operation gestalt, and 
the fluid flow which passes an adjustable converging section can be computed to accuracy. 
[0063] 

[Effect of the Invention] Since the fluid flow is computed without being based on the relational expression 
of adiabatic change according to this invention as explained above, calculation of the fluid flow can 
compute to accuracy under the conditions on which entropy conditions are not satisfied. 
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